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ABSTRACT 

We model the spectral energy distribution of NGC 4258 with updated broadband observa- 
tions from radio to X-rays and a coupled accretion-jet model that surrounding a Kerr black 
hole (BH). The observed radio jet and the warped H2O maser disk are assumed to be trig- 
gered by a spinning BH through Blandford-Znajek mechanism and Bardeen-Petterson effect 
respectively. The accretion flow is modeled as an inner radiatively inefficient accretion flow 
(RIAF) and an outer truncated standard thin disk, where the transition radius is ~ 10 3 i? g 
based on the width and variability of narrow Fe Ka line. We find that the RIAF surrounding a 
mildly spinning BH with dimensionless spin parameter a* ~ 0.73 can well reproduce the ob- 
served X-ray emission and jet power, where the observational jet power is estimated from its 
low-frequency radio emission. The outer thin disk mainly radiates at near infrared waveband, 
while the jet contributes predominantly at radio waveband. We propose that the accretion flow 
initially may consists a two-phase disk that the cold maser disk is sandwiched by substantial 
hot plasma, since that Bondi radius for the diffused hot plasma is much larger than the size of 
the observed maser disk. This can naturally explain why the accretion rate of inner RIAF is 
larger than that of outer maser disk. Using above estimated BH spin and the inferred accretion 
rate in maser disk based on the physical existence of the H2O maser, we find that the warp 
radius of the maser disk that driven by the Bardeen-Petterson effect is ~ 8.8 x 10 4 i? g , which 
is consistent with the observational result very well. 

Key words: accretion,accretion disks-black hole physics-galaxies:nuclei-galaxies:jets- 
galaxies:individual:NGC 4258 



1 INTRODUCTION 

NGC 4258 (M106) is a low-luminosity, weakly active Seyfert 1.9 
(Ho et al. 1997), located at a distance of 7.2 Mpc (Herrnstein et 
al. 1999), where the active nucleus has been detected through radio 
(e.g, Turner & Ho 1994), infrared (e.g, Chary et al. 2000) and X- 
ray observations (e.g, Fruscione et al. 2005; Reynolds et al. 2009). 
Wilkes et al. (1995) also detected the faint blue continuum emis- 
sion and relatively broad (~ 1000 km s , FWHM) emission lines 
in polarized light, which further support the existence of a weakly 
active nucleus. In the radio waveband, high-resolution interfero- 
metric images show the presence of a compact jet (e.g., Herrnstein 
et al. 1997), and two kpc-scale radio spots located at south and 
north of central engine respectively (e.g., Cecil et al. 2000). Fur- 
thermore, VLBI observations have found a set of H2O masers that 
trace a nearly edge-on, geometrically thin gas disk in the region of 
~ 0.14-0.28 pc from the nucleus, where the projected position of 
the maser spots indicate that the maser disk is substantially warped 
(Miyoshi et al. 1995). The nearly perfect Keplerian rotation curve 
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of these H2O masers provides the precise measurement of the black 
hole (BH) mass 3.9 ±0.1 x 10 7 Mq and the inclination angle of 
~ 83° for the outer maser disk (Miyoshi et al. 1995; Herrnstein et 
al. 1999). 

The bolometric luminosity of NGC 4258 is small compared 
its Eddington luminosity (Lboi ~ 10 _4 I/Edd, e.g, Reynolds et al. 
2009, where L Edd = 1.3 x 1O 38 (M BH /M )erg s" 1 ), which in- 
dicates that it is a typical low-luminosity active galactic nucleus 
(LLAGN). Its low luminosity may be caused by the radiative in- 
efficiency of the accretion flow at low accretion rate, where the 
hot, optically thin, geometrically thick radiatively inefficient ac- 
cretion flow model has been developed in last several decades 
(RIAF, an updated version of the original advection dominated ac- 
cretion flow model, ADAF; e.g., Ichimaru 1977; Naraya & Yi 1994, 
1995; see Narayan & McClintock 2008 for a recent review). The 
RIAF model can successfully explain most observational features 
of nearby LLAGNs (e.g., Quataert et al. 1999; Cao & Rawlings 
2004; Nemmen et al. 2006; Wu et al. 2007; Yuan et al. 2009; Xu & 
Cao 2009; see Ho 2008 for a recent review). For accretion rate less 
than a critical value (M cr it — O.OlMEdd), the standard thin disk 
(Shakura & Sunyaev 1973) may transit to a RIAF at inner region. 
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It was found that the outer thin disk with a transition radius i?tr ~ 
several ten to several hundred gravitational radii normally can well 
explain the "red bump" that observed in the spectral energy distri- 
bution (SED) of LLAGNs (e.g., Quataert et al. 1999; Ho 2008; Xu 
& Cao 2009). The spectrum of NGC 4258 has been modeled with 
the truncated disk model, where the X-ray spectrum is explained 
by a RIAF and the infrared spectrum is fitted by the outer thin disk 
(Lasota et al. 1996; Gammie et al. 1999). Based on the possible 
power-law distribution of the infrared emission after correcting 18 
mag of visual extinction, Yuan et al. (2002) proposed that the spec- 
trum from infrared to X-rays can also be explained by a pure jet 
model. 

Although the warped disk of NGC 4258 has been observed 
nearly twenty years, the nature of physical mechanism responsible 
for producing such a structure is still unclear. There are several pos- 
sible explanations for origin of the warped disk. Papaloizou et al. 
(1998) proposed that it could be produced by a binary companion 
orbiting outside the maser disk, where the companion object should 
have a mass comparable to that of disk. However, there is still no 
observational evidence for such an object. The second possibility is 
that the warped disk is caused by the non-axisymmetric forces due 
to the radiation pressure (Pringle 1996). However, this mechanism 
cannot produce the inferred warp in NGC 4258 if the accretion disk 
is radiatively inefficient (e.g., Maloney et al. 1996). For the possi- 
ble radiatively efficient standard disk, it will be stable against the 
radiation instability and will not warp if the viscosity parameter of 
the maser disk a < 0.2 (Caproni et al. 2007). Caproni et al. (2007) 
concluded that the warping of the maser disk in NGC 4258 should 
be caused by the Bardeen-Petterson effect, which is predicted in the 
framework of general relativity when the spin of a Kerr BH is in- 
clined in relation to the angular momentum vector of the accretion 
disk (see also Martin 2008). The Lense-Thirring precession drives 
a warp in the disk, and the inner part of the disk is aligned with the 
BH (Bardeen & Petterson 1975). 

The jet formation mechanisms have been extensively stud- 
ied in the last several decades. The recent magnetohydrodynamic 
(MHD) simulations and analytical works showed that both BH spin 
and large-scale magnetic field may play important roles in jet for- 
mation (e.g., McKinney & Gammie 2004; Hirose et al. 2004; De 
Villiers et al. 2005; Hawley & Krolik 2006; Hawley & Krolik 2006; 
Li & Cao 2009; Tchekhovskoy et al. 201 1; Cao 2012; Cao & Spruit 
2013). In the Blandford-Znajek (BZ) process, energy and angu- 
lar momentum are extracted from a spinning black hole (BH) and 
transferred to a remote astrophysical load by open magnetic field 
lines. Narayan & McClintock (2012) found that the scaling of the 
jet power and the BH spin in X-ray binaries is roughly consistent 
with that predicted by Blandford & Znajek (1977), where the BH 
spins are determined via the continuum-fitting method (e.g., Zhang 
et al. 1997; Gou et al. 2011). Wu et al. (2011) also proposed that 
the radio FR I galaxies are powered by the fast rotating BHs based 
on the BZ process. 

NGC 4258 is the only AGN for which the BH mass, distance 
and the geometry of the outer disk are so accurately known. There- 
fore, the central engine of NGC 4258 have been investigated care- 
fully in last decade. However, both the multiwavelength observa- 
tions of NGC 4258 and theories of the accretion flow have made 
important progress in last several years. For example, the X-ray ob- 
servations of the narrow Fe Ka line give a strong constraint on the 
possible transition of the outer SSD to an inner RIAF at ~ 10 R s 
(Reynolds et al. 2009); the strong variability of millimeter and X- 
ray emission suggest that they should originate in the region near 
the BH (Reynolds et al. 2009; Doi et al. 2013); the mid-infrared 



emission of NGC 4258 is found to similar to those of conventional 
bright Seyferts with no evident jets (Mason et al. 2012), etc. On 
the theoretical side, the global, time-dependent, numerical simu- 
lations revealed that the outflows may be important in RIAFs and 
only a fraction of the mass that available at large radii actually ac- 
cretes onto the BH (e.g., Blandford & Begelman 1999; Narayan et 
al. 2012; Yuan et al. 2012). Furthermore, in the old ADAF model, 
the gravitational energy is assumed to predominately heat the ions. 
However, it was later realized that the electrons may also receive 
an important fraction of the dissipated energy through the magnetic 
reconnection and/or the viscosity arising from anisotropic pressure 
(e.g., Sharma et al. 2007). In this work, we re-explore the central 
engine of NGC 4258 with the updated observations, and we use the 
accretion-jet model surrounding a Kerr BH that is not considered in 
former works. Furthermore, we also investigate whether or not the 
Bardeen-Petterson effect that driven by a spinning BH can produce 
the warped maser disk as observed in NGC 4258. 

For a distance of 7.2 Mpc and a BH mass of 3.9 x 10 7 Mq 
for NGC 4258, 1" corresponds to 35 pc and 1 pc corresponds to 
5.4 x 10 5 i? g , where the gravitational radius R g = GMbh/c 2 . 



2 THE NUCLEAR SED 

The radio emission of NGC 4258 have been thoroughly investi- 
gated from ~ 1 GHz upto millimeter waveband (see Doi et al. 2013, 
and references therein). We simply summarize its radio properties 
as follows. VLBI imaging at ~ 1.5 GHz showed a two-sided nu- 
clear jet extending up to ~ 10 6 -R g (Cecil et al. 2000), which may 
physically connect with the kpc-scale jet-like structures (anoma- 
lous arms, e.g., van der Kruit et al. 1972). Herrnstein et al. (1998) 
reported a jet-like structure in 22 GHz radio observation, which 
oriented along the rotation axis of the maser disk, where the north- 
ern and southern components have flux density of 2.5-3.5 mjy and 
0.5 mjy respectively. Interestingly, the VLBI images also provided 
a 3<t upper limit of 220/^Jy on 22 GHz at the dynamical center, 
which place a strong constraint on accretion disk (Herrnstein et al. 
1998), where the actual upper limit is likely 1-2 times higher than 
the observed 220^iJy if considering the possible free-free absorp- 
tion in the disk (e.g., Herrnstein et al. 1996). Other multi-frequency 
radio observations (e.g., 5, 8.4, 15 and 43 GHz) at ~ 0.5 resolu- 
tion toward the nucleus are also presented in Nagar et al. (2001) and 
Doi et al. (2013). Doi et al. (2013) also presented the millimeter and 
submillimeter observations for NGC 4258, where the strong vari- 
ability at 100 GHz observed by the Nobeyama Millimeter Array 
(NMA) suggest it should come from the region close the BH, while 
most of the submillimeter emission at 347 GHz observed by lames 
Clerk Maxwell Telescope (ICMT) may originate in interstellar dust 
of the host galaxy. In building the radio spectrum, the average flux 
is used if there exist multi-epoch data at given waveband (see Table 
1 and Doi et al. 2013, for details). 

Chary et al. (2000) presented high-resolution infrared imaging 
of the nucleus in NGC 4258 from 1 to 18 /im. The compact source 
of radiation is still unresolved even at near-infrared resolution of 
0.2 . Chary et al. (2000) found that the near- and mid-infrared data 
can be fitted by a single power law as V ~ 1A±0A if the infrared 
data are corrected for about 18 mag of visual extinction, and they 
suggest that the power-law infrared emission may have nonthermal 
origin. However, Mason et al. (2012) found that the mid-infrared 
emission of NGC 4258 is similar to those of conventional bright 
Seyferts, and the quite evident silicate emission features demon- 
strate that the dust should present in its nuclear region. Therefore, 



© RAS, MNRAS 000,[Tp?? 



Warping of accretion disk and launching of jet by a spinning black hole in NGC 4258 3 



Table 1. Spectral energy distribution for NGC 4258. 



Waveband 


Hz 


Flux 
mjy 


Luminosity 

erg s 1 


Resolution 


References 


1.4 GHz 


1.4xl0 9 


2.73 


35.37 


~ l" 


1 


5 GHz 


5xl0 9 


2.0±0.1 


35 7 Q+0.02 
JJ - "-0.02 


~ 0.5" 


2 


8.4 GHz 


8.4xl0 9 


2.2±0.1 


36 06+ 02 


~ 0.5" 


2 


15 GHz 


1.5 xlO 10 


2.7±0.2 


36 40+0.03 


~ 0.5" 


2 


22 GHz 


2.2 xlO 10 


2.9±0.4 


36 65+ ' 05 


~ 0.5" 


2 


22 GHz 


2.2 xlO 10 


<0.22 


< 35.48+ 3 


~ 0.0003" 


3 


43 GHz 


4.3 xlO 10 


6.9±0.2 


37 26+ - 01 
J '- ZD -o.oi 


~ O.l" 


2 


100 GHz 


lxlO 11 


6.9±0.2 


37 63+ ' 01 


~ 10" 


2 


347 GHz 


3.47x10" 


93.7±20.8 


39.301°?? 


~ 20" 


2 


17.9/xm 


1.68xl0 13 


300±30 


41 49+0-05 
^ 1>fy -0.06 


~ 0.5" 


4 


12.5/xm 


2.40xl0 13 


165±20 


41 S 9+0.06 


~ 0.5" 


4 


10.5/im 


2.86xl0 13 


100±12 


41.25+°' i 

—0.06 


~ 6" 


4 


3.45/xm 


8.70xl0 13 


20 ±3 


41-03+n? 


~ 0.2" 


4 


2.2/^m 


1.36xl0 14 


4.0±0.7 


40.53±C;?1 


~ 0.2" 


4 


1.6/^m 


1.88xl0 14 


0.9±0.3 


40 n?+° 12 


~ 0.2" 


4 


1.1 (im 


2.73xl0 14 


<0.5 


<39.93 


~ 0.2" 


4 


B-band 


6.81xl0 14 




39.97 




5 


2-10 keV 






40 Q4+ 016 




6 



References: (1) Ho & Ulvestad (2001); (2) Doi et al. (2013); (3) Herrnstein et al. (1998); (4) Chary et al. (2000); (5) this work; (6) Frascione et al. 
(2005). 



the possible power-law distribution of infrared emission as pro- 
posed by Chary et al. (2000) may be not intrinsic. The observed 
high-resolution infrared data without any corrections are listed in 
Table 1 (see Chary et al. 2000 for more details). 

The optical emission of NGC 4258 is less clear, due to the 
possible obscuration by the putative torus or/and warped disk. The 
indirect estimation of the optical luminosity is L ° ~ 10 37 — 

r J 5500A 

10 39 erg s _1 A _1 from the nuclear polarized light (Wilkes et al. 
1995). However, this estimation still suffer many uncertainties. In 
this work, we estimate the nuclear optical emission of NGC 4258 
from its H/3-line luminosity, using the correlation between in/3 and 
Mb for both bright and faint type I AGNs, as calibrated by Ho & 
Peng (2001). From this method, we find that its B-band luminosity 
is Lb = 10 39 ' 97 erg s _1 , where the distance-corrected H/3-line 
luminosity Lh/3 = 10 38 ' 82 erg s _1 and the empirical correlation 
logL H/3 = (-0.34±0.01)M S + (35.1 ±0.25) were adopted (Ho 
&Peng 2001). 

The X-ray emission of NGC 4258 is known to be variable 
on timescale of hours to years with unabsorbed 2-10 keV X-ray 
luminosity of (4.2-17.4)x 10 40 erg s _1 , where the average X-ray 
luminosity is 8.8 x 10 40 erg s" 1 and the average photon index is 
1.6 ± 0.2 (see Frascione et al. 2005). Using the data of Suzaku, 
Reynolds et al. (2009) proposed that the region of the hard X-ray 
emission should be smaller than 257? g based on the high-amplitude 
intraday variability. They further constrained that the narrow Fe Ka 
line emitting region to be ranged from ~ 3 x 10 3 i? g to ~ 4 x 10 4 i? g 
based on the line width and the flux variability by comparing the 
Suzaku data with XMM-Newton data that taken 160 days later. 

All the data for building the SED of NGC 4258 and related 
references are listed in Table 1. 



3 MODEL FOR THE CENTRAL ENGINE 

Our model consists three components: an inner RIAF, an outer trun- 
cated thin disk that may connect to the warped maser disk and a jet, 
which are illustrated in Figure 1. We describe these components in 
details as follows. 



3.1 The inner RIAF 

The RIAF model surrounding a Kerr BH is considered in this work. 
Only main features of the model are described here, and the de- 
tails of which can be found in Wu & Cao (2008, and references 
therein). We employ the approach suggested by Manmoto (2000) 
to calculate the global structure of the RIAF in the general rela- 
tivistic frame, which allows us to calculate the structure of a RIAF 
surrounding either a spinning or a nonspinning BH. Both observa- 
tions and theories suggest the existence of outflow in hot accretion 
flows. To describe the possible mass loss of the disk, we assume 
that the accretion rate is a function of radius, M = M t i(R/Rti) a , 
where R tr is the outer radius of the RIAF and M tr is the accre- 
tion rate at R tT (e.g., Blandford & Begelman 1999). Yuan et al. 
(2012) found that the profile of M(R) has two parts from the sim- 
ulations. When R < 10R g , s ~ 0; but outside 10i? g , s ~ 0.5. 
In the present paper, we set s = 0.4 throughout the flow for sim- 
plicity, which is also roughly consistent with that constrained from 
the observation of the supermassive BH in our Galactic center-Sgr 
A* (Yuan et al. 2003). The global structure of the RIAF surround- 
ing a BH that spinning at a rate a* with mass Mbh can be cal- 
culated with proper outer boundaries, if the parameters m, a, ft 
and S are specified (e.g., Manmoto 2000), where all radiation pro- 
cesses (Synchrotron, Bremsstrahlung and Compton scattering) are 
included self consistently. The parameter a* = J/GM^c -1 is 
the dimensionless angular momentum of the BH (J is the angu- 
lar momentum of the BH), and rh = M /A^Edd is the dimen- 
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R w ~10 5 R a 




R B ~10 6 R g 

R n .~10 3 R s 



Figure 1. The cartoon illustrating the model for the central engine of NGC 
4258. It mainly consists of three components: an inner RIAF, an outer trun- 
cated thin disk which may connect to outer warped maser disk, and a jet. 
The dotted lines represent the Bondi radius, and the hot gas inside this ra- 
dius will also be captured by BH. 

sionless accretion rate (the Eddington accretion rate is denned as 
M E dd = 1.3 x 10 18 Mbh/Mq g s" 1 ). For the viscosity parame- 
ter a, we adopt the typical value of a — 0.2 which is constrained 
from the dwarf novae in hot state or the modeling the spectrum of 
LLAGNs (a ~ 0.1 - 0.3, Narayan & McClintock 2008, and refer- 
ences therein). Numerical simulations suggest that magnetic fields 
are generally subthermal, with P mag /P gas ~ 0.1 (Hawley et al. 
2006), so we expect the ratio of the gas to the total pressure (sum of 
gas and magnetic pressure) ft = P g /Ptot = 0.9. The parameter 8 
describes the fraction of the turbulent dissipation that directly heats 
the electrons in the flow. In this work, we adopt S = 0.1, which can 
well describe the multiwavelength radiation of LLAGNs (Liu & 
Wu 2013). The field-enhancing effect caused by frame dragging is 
also considered (Meier 2001). The amplified magnetic field related 
to the magnetic field produced by the dynamo process in the RIAF 
can be expressed as B = gBdynamo, where g — Q/fl is the field 
enhancing factor, the disk angular velocity fHs a sum of its angular 
velocity relative to the local metric f2 plus the angular velocity of 
the metric itself in the Boyer-Lindquist frame u> = —g^tjg^, i.e., 
fl = Q. + u). The spectrum is calculated using the global structure 
of RIAF surrounding a Kerr BH, where the gravitational redshift 
effect is considered, while the relativistic optics near the BH is ne- 
glected (e.g., Manmoto 2000). 

3.2 The outer thin disk 

Outside the inner RIAF there is an outer thin disk with a truncation 
radius of R tr . The physical mechanism for the disk transition is 
still an open issue. In this work, we set R tr = 1000P g for NGC 
4258, which is constrained from the width and flux variability of the 
narrow Fe Ka line (Reynolds et al. 2009). The outer radius of the 



thin disk is unclear, which may smoothly connect to the observed 
maser disk. 

In addition to the accreting cold gas through the maser disk 
at large scale, the hot gas around the Bondi radius also provide a 
suitable mass supply rate for BH accretion (e.g., M 87, Di Matteo 
et al. 2003). The Bondi radius R B = GM B h/c 2 ~ 10 6 P g for 
a typical temperature of 1 keV for the diffused hot gas at ~ Rb, 
where c s ~ 0.1T 1 ^ 2 km s _1 is the sound speed of the hot gas. 
The Bondi radius is around one order of magnitude larger than the 
observed outer radius of maser disk (~ 0.28pc ~ 1.5 x 10 5 R g ). 
Therefore, the inflow initially should consist a two-phase accretion 
flow that the cold thin disk is sandwiched by hot plasma. The hot 
inflow inside the Bondi radius will be gradually cooled down and 
condensed into the thin disk due to the increase of the radiative ef- 
ficiency through the inverse Compton process when the cold disk 
is present. Based on this scenario, we think that it is possible for 
the accretion rate at truncation radius, A/ tr , larger than that of the 
maser disk, M m . In this work, we constrain Mtr through modeling 
of the X-ray spectrum and jet power of NGC 4258, and then check 
the radiative spectrum of the outer SSD with Mtr and R tI . Most 
of the radiation from the truncated thin disk originate in the region 
near the transition radius due to the 1/72 dependence of the gravita- 
tional potential. Therefore, we calculate the spectrum of outer thin 
disk simply using a constant rate at R tT , which will not affect our 
main conclusion. The thin disk emits locally as a blackbody, and its 
spectrum can be calculated if BH mass, accretion rate and inclina- 
tion angle of the disk are known. We adopt the inclination angle of 
the inner disk to be 74° that may aligned with the equatorial plane 
of the BH (e.g., Krause et al. 2007). 



3.3 The jet: power & spectrum 

We investigate both the total power and radiative spectrum of the jet 
in NGC 4258, since that the radio jet has been detected directly. To 
evaluate the jet power extracted from the inner region of RIAF that 
surrounding a spinning BH, we adopt the hybrid model as a variant 
of a BZ model proposed by Meier (2001). The total jet power is 
given by 

Pjet = B*R 4 n 2 /32c, (1) 

where R is the characteristic size of jet formation region, the 
poloidal magnetic field B p ~ <?Pdynamo (see Section 3.1 and 
Meier 2001 for more details). Following the work by Nemmen et al. 
(2007), all the quantities are evaluated at the innermost marginally 
stable orbit of the disk R = P ms , which is roughly consistent with 
the jet launching region in the numerical simulations. Wu et al. 
(201 1) found that the jet efficiency of this model is consistent with 
those of MHD simulations very well. 

Both the protons and electrons will be accelerated by the 
magnetic field energy that extracted from the disk and/or rotating 
BH, and the kinetic energy of particles will be eventually released 
through the possible shock in the jet. To describe the jet radia- 
tion, we adopt the internal shock scenario that has been used to 
explain the broadband SED of X-ray binaries, AGNs and afterglow 
of gamma-ray burst (e.g., Piran 1999; Yuan et al. 2005; Wu et al. 
2007; Nemmen et al. 2012). The jet is assumed to have a conical ge- 
ometry with opening angle <j>. The internal shock in the jet occur as 
a result of the collision of shells with different velocities, and these 
shocks accelerate a small fraction £ c of the electrons into a power- 
law energy distribution with index p. From the shock jump condi- 
tions, the electron and ion temperatures can be derived. Following 
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the widely adopted approach in the study of GRBs, the energy den- 
sity of accelerated electrons and the amplified magnetic field are 
described by two free parameters e c and ea , which roughly follow 

1/2 

e c = e B ' (Medvedev 2006). We then calculate the radiative trans- 
fer by both thermal and power-law electrons in the jet, although we 
find that the latter plays a dominant role. Only synchrotron emis- 
sion is considered, since Compton scattering is not important in this 
case (e.g., Markoff et al. 2001; Wu et al. 2007). 

The 22 GHz radio observations show that NGC 4258 have a 
broad jet, where the jet open angle ~ 20° (Herrnstein et al. 1997), 
which is similar to the jet in other nearby LLAGNs (e.g., M81, 
Markoff et al. 2008), and hence we adopt tf> = 20°. The typical 
value of e c =0.1 for LLAGNs is adopted in this work (e.g., Yuan et 
al. 2009; Nemmen et al. 2012), and then 6b=0.01 according to the 
relation proposed by Medvedev (2006). Shock acceleration theory 
typically gives 2 < p < 3, and we fix p = 2.5. For the fraction of 
electrons that were accelerated into the power-law distribution, we 
use £ e = 0.1 as that of FR I radio galaxies (Wu et al. 2007). There 
are two free parameters in the jet model: the outflow rate A4j e t and 
the jet velocity Vj ct , and these two parameters can be determined by 
fitting the spectrum and the jet power that estimated from equation 
(1). 



4 RESULT 

4.1 Estimating the jet power 

We estimate the jet power of NGC 4258 from its low-frequency ra- 
dio emission since it usually originates in the diffused optically thin 
lobes that move at low velocities and do not suffer Dopper boost- 
ing effect. Wu et al. (2011) derived an empirical correlation be- 
tween 151 MHz radio luminosity and jet power that measured from 
X-ray cavities for a sample of low-power radio galaxies, which is 
log P jot = 0.52 ± 0.04 + 22.31 ±1.75 with an intrinsic scatter of 
0.23 dex. Through this correlation and the 151 MHz radio luminos- 
ity of 37.28 erg s" 1 reported in Hales et al. (1988), we find that jet 
power of NGC 4258 is 10 41 ' 69 erg s" 1 . The estimated jet power 
is consistent with that reported in Falcke & Biermann (1999) very 
well, where the jet power of 10 41 ' 7 erg s _1 is derived from the 22 
GHz nuclear radio emission and offset of the core from the dynam- 
ical center. 



4.2 Modeling the multiwavelength spectrum 

The SED of NGC 4258 is plotted in Figure 2. In spectral fitting, 
we, firstly, adjust the parameters of BH spin, a*, and the accretion 
rate at the transition radius, m tr , to fit the observed X-ray spec- 
trum and simultaneously allow the derived jet power from equation 
(1) to be consistent with the above estimated jet power, since that 
both jet power and high-energy X-ray emission are predominantly 
originate from the inner region of the RIAF near the BH. We find 
that the best-fit parameters are a* = 0.73 and m tr = 6 x 10 -3 . 
It is interesting to note that the emitted spectrum from the underly- 
ing RIAF (dashed line) can not only reproduce the observed X-ray 
spectrum fairly well; but also roughly consistent with the estimated 
nuclear optical emission and the upper limit of 22 GHz emission 
that radiated from the central engine (red solid and open circles in 
Figure 2). Secondly, we show the emitted spectrum from the outer 
thin disk in Figure 2 (dot-dashed line), where the inner radius of 
Rtr = 10007? g and the accretion rate rh tr are used directly in spec- 
tral calculations. We find that the outer SSD mainly radiates at near 
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Figure 2. The spectral energy distribution of NGC 4258. The dashed, dot- 
dashed and dotted lines show respectively the emission from the inner 
RIAF, outer SSD and jet, while the solid line shows the sum of the radiation 
from these components. The red solid circle represent the nuclear optical 
emission estimated from the H/3-line luminosity, and the red open circle 
represent the upper limit of the 22 GHz radio emission that coincident with 
the central engine. The blue solid line denote the MIR spectrum of NGC 
4258 that chosen from Mason et al. (2012). 



infrared waveband, and the theoretical prediction is higher than that 
of the observations, where the observed infrared data is not cor- 
rected for the possible obscuration. Thirdly, the radio emission of 
NGC 4258 should originate in the jet since it was observed directly. 
We find that the jet model with outflow rate, mj ct = 5 x 10 -6 and 
jet velocity, Vj c t = 0.86c, can well reproduce both the observed ra- 
dio spectrum (dotted line) and the estimated jet power, where about 
~ 0.83% of the mass that ultimately accreted by the black hole is 
channeled into the jet. The solid line in Figure 2 is the sum of the 
spectrum from the inner ADAF, outer thin disk and the jet. We find 
that our model gives a quite nice fit for the observed spectrum ex- 
cept the mid-infrared emission, which may not originate from the 
central engine of accretion/jet. 



4.3 The Bardeen-Petterson effect and warp radius of outer 
maser disk 

Frame dragging produced by a Kerr BH causes precession of a par- 
ticle if its orbital plane is inclined in relation to the equatorial plane 
of the BH, which is known as Lense-Thirring precession. There ex- 
ist a critical radius, ii w , in the warped disk, where the disk material 
interior to ii w will be aligned with the equatorial plane of the black 
hole, while the outer portion of the disk (R > 7? w ) will maintain its 
original orientation (Bardeen & Petterson 1975). The critical warp 
radius R w can be determined by comparing the angular momentum 
of the spinning BH and the angular momentum of the disk in per 
logarithmic interval of the radius, which is given by Lu & Zhou 
(2005), 



Ryu 
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71 .. in 3 16/35 4/7 
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The motion of the infalling gas of NGC 4258 cannot be measured 
directly, and there are some indirect methods to estimate the ac- 
cretion rate in the maser disk. The maser amplification is favored 
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when the particle density range from 10 8 to 10 10 cm 3 (Moran 
et al. 1995). Based on the standard thin disk model, Cao & Jiang 
(1997) derived an accretion rate of M m ~ 8 x lO~ 4 aM0yr _1 
at the inner edge of the maser disk in NGC 4258 , where they as- 
sume the disappearance of the maser disk at the inner region is 
caused by the density of the molecular gas in disk reach the upper 
limit for the existence of the maser. Using above estimated BH spin 
a, = 0.73 and accretion rate in maser disk, we find that the warp 
radius i? w ~ 1.0 x 10 5 a 2/35 i? g . For the typical value of a ~ 0.1, 
the warp radius R w ~ 8.8 x 10 4 J? g , which is evidently not sensi- 
tive to the viscosity parameter. 



5 DISCUSSION 

The observational features of NGC 4258, such as the hard X-ray 
spectrum, narrow Fe Ka line, infrared bump, and the variable 
inverted radio spectrum, indicate that it is most likely a typical 
LLAGN (Ho 2008). NGC 4258 is also one of the nearest LLAGNs 
that suitable for studying the properties of their central engines, 
following the extremely sub-Eddington accreting source-Sgr A* at 
the Galactic center (8 kpc, Reid 1993) and the closest LLAGN- 
M81 (3.6 Mpc, Freedman et al. 2001). The warped maser disk at 
sub-parsec scales in NGC 4258 places strong constraints on both 
its BH mass, disk geometry and even the possible accretion rate in 
the maser disk. We investigate the properties of the central engine 
in NGC 4258 with a accretion-jet model that surrounding a Kerr 
BH, where we find that a mildly rotating BH with a, ~ 0.73 can 
well reproduce its nuclear SED, jet power and warp radius of the 
maser disk. 

The inner RIAF+outer thin SSD+jet has been considered as 
the standard model for nearby LLAGNs and low/hard state X-ray 
binaries (e.g., Quataert et al. 1999; Nemmen et al. 2006; Xu & Cao 
2009, and see Ho 2008 for a recent review), where the inner RIAF 
mainly radiates at hard X-ray band, the outer SSD radiates at in- 
frared to optical band, while the jet mainly radiates at radio wave- 
band. Our detailed modeling of NGC 4258 do support this scenario. 
It is interesting to note that the RIAF can also well reproduce the 
nuclear optical emission and the upper limit of the compact 22 GHz 
emission in NGC 4258, where the nuclear optical emission is esti- 
mated from the emission-line luminosity and the 22 GHz emission 
originate in the compact central engine. The 6.4 keV Fe Ka line 
(e.g., width, strength, and variability) is one of the most powerful 
probes for the cold gas surrounding the central BH. Assuming the 
Keplerian orbits in the outer SSD, Reynolds et al. (2009) estimated 
the line emitting region roughly range from 3 x 10 3 to 4 x 10 4 i? g 
using the line width and variability, which suggest that the outer 
SSD may transit to an ADAF at R tv < 10 3 R g . Using this tran- 
sition radius and the accretion rate m tr that constrained from the 
fitting of the X-ray emission, we find that the outer SSD mainly 
radiates at near infrared waveband. The theoretical prediction is 
higher than the observed value, which may caused by the possible 
obscuration of the infrared emission by the putative torus or/and 
the warped disk. The exact obscuration is still unclear, which is be- 
yond the scope of this work. The mid-infrared emission (> 10/im) 
is much higher than that predicted by the outer SSD, which sug- 
gest that the mid-infrared emission may come from other com- 
ponents. Laine et al. (2010) found that the mid-infrared emission 
of NGC 4258 at 8/im that obtained from the Spitzer is mainly 
originate from the polycyclic aromatic hydrocarbon molecules and 
hot dust. The well-known MIR— hard X-ray correlation for radio- 
quiet luminous AGNs can also be used to explore the possible ori- 



gin of the MIR emission, where the MIR emission of these radio- 
quiet AGNs should originate from the molecular gas of the putative 
torus since there are no relativistic jets in them. Using the relation 
logLg Mm = 1.14 log I/2-iokcV — 5.40 that constrained from the 
bright Seyferts and quasars (Mason et al. 2012) and the average X- 
ray luminosity of log I/2-iokcV = 40.94, we find that this relation 
predicts log Lsum = 40.27, which is consistent with the observed 
MIR emission, log Lio. 5ftm = 40.25 for NGC 4258 very well. This 
result further support that the MIR emission of NGC 4258 may 
originate in the molecular cloud and/or dust as that of radio-quiet 
luminous AGNs. Chary et al. (2000) proposed a nonthermal origin 
for infrared emission from 1.1 to 18 jj,m based on the power-law 
distribution of the spectrum after correcting for about 18 mag of 
visual extinction. However, the intrinsic extinction is still unknown 
yet and the MIR emission may have a different origin compared the 
NIR emission as discussed above. 

The accretion rate in outer maser disk has been explored with 
different methods, which is ~ (1 — 8) x 10 _4 aMQ/yr, where 
the rate ~ 10 -4 aMQ/yr is estimated by assuming a transition 
from a molecular to an atomic disk at outmost region of maser disk 
(Neufeld & Maloney 1995; Herrnstein et al. 2005); and the rate of 
~ 8 x 10 -4 o:Mq /yr is estimated by assuming the quenched maser 
at its inner region is due to the molecular density reach its upper 
limit of 10 10 cm -3 for the existence of the H2O maser (Cao & 
Jiang 1997). It can be clearly found that the estimated accretion rate 
in maser disk is much lower than that constrained from the inner 
RIAF, where M tr ~ lO~ 2 M /yr (see also Gammie et al. 1999). 
We note that this is possible since that the hot gas inside the Bondi 
radius will also be accreted beside the accreting cold gas through 
the maser disk, where the Bondi radius, Rb ~ 10 6 i? g , is larger 
than the outer radius of the maser disk. The high-resolution X-ray 
observation on the diffused gas near the Bondi radius will help to 
test this issue. The accretion rate of the outer thin disk may be not 
a constant considering the mixture of inflowing hot and cold gas, 
where the hot gas will be condensed down to the SSD gradually due 
to the cooling. Actually, the estimated accretion rate in the inner 
edge of maser disk (Cao & Jiang 1997) is higher than that in outer 
edge (Neufeld & Maloney 1995; Herrnstein et al. 2005), which may 
more or less support this scenario. In calculating the spectrum of 
outer SSD, we simply use the rate at transition radius, which will 
not affect our main conclusion since that most gravitational energy 
of the thin disk is released in its inner region. 

Both the radiation of the RIAF and the jet power extracted 
from the disk/BH increase with the increasing of the BH spin, 
where, however, jet power is much more sensitive to the spin pa- 
rameter compared the X-ray spectrum. To reproduce the observed 
X-ray spectrum and the constrained jet power simultaneously, we 
find that the BH should be mildly spinning with a* — 0.73 in the 
NGC 4258. We find that the X-ray emission from the RIAF is sen- 
sitive to the parameter of 8 and less sensitive to other parameters. 
The jet power that extracted from the disk/BH is less affected by 
model parameters for given BH spin and the accretion rate (Wu & 
Cao 2008, and reference therein). Sharma et al. (2007) found that 
the parameter 8 may be in range of ~ 0.03 — 0.3 based on the 
simulations, depending on the model details. To reproduce the ob- 
served X-ray luminosity and the estimated jet power, we find that 
the spin should be a, ~ 0.85 for 8 — 0.3, where the accretion rate 
of the inner RIAF is slightly lower due to more dissipated energy 
is directly heated electrons and, therefore, the slightly higher BH 
spin is needed to reproduce the estimated jet power. For the case 
of lower value of 8 = 0.03, we find the BH spin should be also 
slightly lower with a* ~ 0.6. It can be found that the BH is mildly 
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spinning for both cases even considering the possible uncertainties 
of the parameter 8. However, we find that the RIAF spectrum with 
8 ~ 0.1 cannot only reproduce the observed X-ray spectrum; but 
also roughly consistent with the estimated optical emission and ob- 
served 22GHz compact radio emission, where all these emission 
may originate from RIAF. Therefore, a, ~ 0.73 may be the best- 
fit parameter in our model. We adopt the shock scenario of the jet 
model to fit the radio emission of NGC 4258, since the radio is 
emitted predominantly by the jet. In spectral fitting, we find that 
the ratio of the mass loss rate in the jet to the accretion rate near the 
horizon of the BH is ~ 1%, which is a typical ratio for BH that ac- 
creting at a low rate (Yuan & Cui 2005; Yu et al. 201 1; Nemmen et 
al. 2012). Using this outflow rate, we find that the speed of the jet 
is ~ 0.86c, which is mildly relativistic. Considering the possible 
slightly lower fraction of electrons are accelerated into the power- 
law distribution, e.g., £ c = 0.01, in less powerful jet of NGC 4258 
compared that of £ e = 0.1 in the powerful jets of FR Is, we find that 
the mass outflow rate will be slightly higher, mj ct = 2.0 x 10~ , 
and the outflow speed will be lower, Vj ct = 0.6c, for given esti- 
mated jet power. Therefore, the jet in NGC 4258 should be mildly 
relativistic even considering possible uncertainties in model param- 
eters. In fact, the mildly relativistic jets have been inferred from 
observations of other LLAGNs/nearby Seyfert galaxies, e.g., Arp 
102B (0.45c; Fathi et al. 2011) and NGC 4278 (0.76c; Giroletti 
et al. 2005). Most importantly, we find that the jet model cannot 
reproduce the radio and infrared/X-ray emission simultaneously, 
even we allow all parameters of jet model to be free. Therefore, 
the nuclear SED of NGC 4258 is most probably produced by the 
coupled accretion and jet. 

The warp of outer maser disk could has various causes. 
Perhaps, the most promising mechanism should be the Bardeen- 
Petterson effect that caused by the interaction of the accretion disk 
and the spinning BH if their spin axes are misaligned (see Intro- 
duction). The superposed maps of maser emission in NGC 4258 
for all epochs suggest that the warp radius should be around 4 mas 
(~ 8 x 10 4 i? g ), where the innermost maser spots at ~ 3 mas deviate 
from the others evidently (see Figure 6 in Argon et al. 2007). We 
find that the critical radius of the warp disk, i? w ~ 8.8 x 10 4 i? g , 
estimated from the Bardeen-Petterson effect (equation 2) is consis- 
tent with the observational results remarkably well, where the BH 
spin parameter is estimated in this work and the accretion rate of 
the outer maser disk is inferred from the physical condition for the 
existence of the maser disk (Cao & Jiang 1997). This result support 
that the Bardeen-Petterson effect should be the physical reason for 
the warped disk in NGC 4258. 
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